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Abstract 


The hyperbolic partial differential equation (PDE) has important practical 
uses in science and engineering. This article provides an estimate for solving 
the Goursat problem in hyperbolic linear PDEs with variable coefficients. 
The Goursat PDE is transformed into a second kind of linear Volterra in- 
tegral equation. A convergent algorithm that employs Taylor polynomials 
is created to generate a collocation solution, and the error using the maxi- 
mum norm is estimated. The paper includes numerical examples to prove 


the method’s effectiveness and precision. 
AMS subject classifications (2020): 45D05, 34K28, 45L05, 30K05. 


Keywords: Hyperbolic partial differential equations; Goursat problem; Volterra 


integral equation; Collocation method; Taylor polynomials. 


1 Introduction 


Partial differential equations (PDEs) are commonly used to address problems 
in various fields, such as engineering, physics, and finance. These equations 
are crucial in the study of various phenomena such as electric currents, grav- 
ity, heat transfer, water wave movement, fluid mechanics, electromagnetism, 
elasticity, quantum mechanics, population dynamics, stock and option pric- 
ing, chemical reaction-diffusion, as well as in the modeling of Schrédinger’s 
equation. 

A second-order linear hyperbolic PDE with two variables for w(e,7) has 


the form 


O?w 02w O2w Ow Ow 
Axa + Bag, + Cop t Dae + Bg, + Pet =0, (1) 


where B? — 4AC > 0, and A,B,C,D,E,F, and G are functions of the 
variables « and 7. By a suitable change of the independent variables, we 
shall show that any equation of the form (1) can be reduced to the canonical 
form or normal form of the hyperbolic equation called the Goursat problem 


oPalt,s) (a dw(t, s) att) 


atas at’ as (2) 
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Equations (1) and (2) are hyperbolic PDEs and appear frequently in the study 
of ecological and cosmological phenomena [26]. To accurately simulate these 
events, it is important to develop efficient and reliable methods for solving 
these equations. Several numerical techniques have been proposed for this 
purpose, including finite difference methods [14, 4], finite element methods 
[2], Taylor matrix method [7, 6], Legendre multi-wavelet Galerkin method [28, 
27], rational Chebyshev method [10], modified variational iteration method 
[1], and Chebyshev wavelet scheme [13]. 


The numerical solution of Goursat problem (2) has been investigated by 
many authors. For instance, Scott [25] considered the homogeneous Goursat 
problem (3) with the coefficients depending on the same variable. Evans 
and Sanugi [12] proposed a nonlinear trapezoidal formula based on geometric 
means. Day [9] used the Runge-Kutta method to approximate the solution 
of (2). Pandey [24] presented a novel exponential finite difference to obtain 
a numerical solution of (2). Drignei [11] developed an algorithm to find the 


quadruple solution to a Goursat problem in a triangular domain. 


In this paper, we consider a collocation solution of the linear Goursat 


problem of the second-order with variable coefficients 


u(t, 8) dw(t, 8) Aw(t, 8) 7 
atOs ae + Olt) SF elt, s)w(ts 8) = f(t,8), (3) 


w(0, s) = a(s), w(t, 0) = B(t), a(0) = B(0), (t, s) € [0, X] x (0, Y], 


+ a(t, s) 


where a, b, c, and f are smooth functions through the domain of discussion. 


For the existence and uniqueness of the solution, see [15]. 


The Goursat PDE (3) is converted to the linear Volterra integral equation 
(VIE) of the second kind 


z 


w(T, z) =9(T,2) + [ k(t, z)w(t, z)dt +f K2(T, 8)w(T, s)ds 


0 


“ | , i ' - sus Basak, (4) 


where g, 1, K2, and «3 are defined in (7). Moreover, the special case of (4) for 
k1 = K2 = 0 is considered in [19]. We develop the collocation method intro- 
duced in [3, 16, 17, 18] to solve (4). This method is based on approximating 


the exact solution of a given integral equation with a suitable function be- 
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longing to a chosen finite-dimensional space (8) such that the approximated 
solution satisfies the integral equation on collocation points (9). An impor- 
tant aspect of this method is that the number of subintervals and the degree 
of Taylor polynomials can be changed to get the best possible result. It is 
also easy to implement, and the approximate solution is based on iterative 
formulas without needing to solve any algebraic equations. 

The paper is organized as follows: In section 2, the converting of PDE 
to the two-dimensional VIE. In Section 3, the approximating solution of (6) 
in each collocation point by a Taylor polynomial. The convergence analysis 
is investigated in section 4. Numerical examples are provided in section 5 to 


illustrate the theoretical results. Finally, a conclusion is given in section 6. 


2 Converting PDE to two-dimensional VIE 


In this section, we study the technique that will convert PDE (3) to an 
equivalent VIE (4). 
Integrating both sides of (3) with respect to s yields 


Ow(t,z)  Ow(t,0) i Fe) eae 2 i: H(t, 9 eas 
0 0 2 


ot ot ot 


+f o(t.s)u(t,s)ds= | f(t, s)ds 
0 0 
which implies 


ff ste.s)as =) ett.0) + [ate Ot Sas La ues 
ab 


~ b(t, 0)w(t, 0) — a ois) re sacfaenen 


(5) 


Integrating again both sides of (5) with respect to t yields 


[ [ f(t, s)dsdt =w(r, z) — w(0, z) — w(7, 0) + w(0, 0) 


+f fa. 3) eS) sar 


A a b(t, z)w(t, z)dt — [ b(t, 0)w(t, 0) dt 
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- ff Aue sjasdt + [fo (t, s)wu(t, s)dsdt. 


[ . f(t, s)dsdt =w(r, z) — a(z) — B(r) + w(0, 0) 
+f a(T, 8)w(T, sjas+ | b(t, z)w(t, z)dt 
-| a(0, s)a(s a b(t, 0) 8 (E)dt 
0 


-[ [= w(t, s)dsdt 
- ff ee (raddsar+ ff ol (t, s)w(t, s)dsdt, 


which is equivalent to the following 2D-VIE: 


Hence, 


z 


w(t, z) =9(T,2) + a Kx (t, z)w(t, z)dt +f k2(T, 8)w(T, s)ds 


0 
+f i k3(t, s)w(t, s)dsdt, (6) 
where 
ki(t,z): = —D(t, z), 
ko(T, 8) : = —a(rT, 8), 
k3(t,s):= alt 2 + ae s) c(t, s), 
g(t, 2) : = a(z) + B(r) — w(0,0) + [ a(0, s)a(s)ds 
+f b(t, 0) 8 (t)dt +f | f(t, s)dsdt. (7) 


3 Taylor collocation method 


In this section, we approximate solutions of 2D-VIE (6) in the space 
Si) ly, a) = {0 : Pam = AUDa im € Mp1 = 0,1,...,N-1m=0,1,...,M—1} (8) 


of the real bivariate polynomial spline functions of degree (at most) p—1 in 
7 and z. Its dimension is NMp?. Here, IIy = {7; = ih,i = 0,1,..., N} and 


Tran. J. Numer. Anal. Optim., Vol. 14, No. 2, 2024, pp 613-637 


Birem, Boulmerka, Laib and Hennous 618 


II = {z; = jk,j =0,1,..., MM} denote, respectively, uniform partitions of 
the intervals [0,7] and [0, Z] with the step-sizes given by h = 4 and k = 4. 
These partitions are define a grid for D 


lina = Uy * Dye = [ta 2m), 0 Sn < NO m <M}. (9) 


Set the subintervals 
On = [T3 Tn41)s7 = 0, 1, Sere /N = 2; ON-1 = [Ty-1, TN]; 


Om = [2mj%m41),m =0,1,...,M—2;  dv_i = [zm-1, 2m], 
and Dyim i= On X om for alln =0,1,...,N—-—1;m=0,1,...,M-—-1. 


To define the collocation solution, we use the Taylor polynomial on each 
rectangle Dpm;n = 0,1,...,N —1;m = 0,1,...,Mé —1. Note that the 
solution w of (6) is known at point (0,0): w(0,0) = g(0,0). 


First, we approximate w in the rectangle Doo by the polynomial 


e 1 d'tiw(0,0) 


v0,0(T; 2) aa 2» a oe ° = ; (7, Zz) E Doo; (10) 
i+j=0 
tg i+j 
where oel0.0) is the exact value of —_ at point (0,0). 


Differentiating equation (6) j-times with respect to z and i-times with 
respect to 7, we obtain 


it g Anas 
0 , one) =9 a) g(r, 2) 
T*OZ 
j i-1 2 : ; 
G\ (i — Vy OPT Fy att la(r, z) 
+ eg (| ae 
XG) n Vee cS al Or" Oz 
j-1 i ; ; ; 
r\ fi, tN [ Qi-1—l antta(r, 2) 
oe (1) () Ori-n aa 2| arnazl 
1=0 n=0 
gj-li-1.. : . ‘ 
g—lyfi-1ly 0-1-7 [ as-1-! Ontlas(z, z) 
rT ( l )( n == pao) peas 
1=0 n=0 
Second, we approximate w in the rectangles D,o,n = 1,...,N —1 by the 
polynomials 


p-1 sii 
1 OD 5 ola) i.j 
Un,o(T, 2) = S- il ridzi (tT -T™)'2z? 3 (7,2)€ Dao, (12) 
i+j=0 
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where @,,,9 is the exact solution of the integral equation: 


z 
On,o(T; z) =9\T, z) +/ K2(T, 8)on,0(T, s)ds 
0 
mel testi T 
+> / rilt2veo(tz)at + f K1(t, z)On,0(t, z)dt (13) 
é=0°7E Tm 
n—1 TE41 z iF z 
a > / / K3(t, 8)v¢,0(t, s)dsdt +/ K3(t, 8)On,0(t, s)dsdt, 
€=0"7E 0 ™m V0 
and 2'**®n.o(™9) is the exact value of 2"7%"0 at point (Tm; 0) 
Or*OzI Or*OzI Pp mp 


Differentiating equation (13) j-times with respect to z and i-times with 
respect to 7, we obtain 


O°+9 bn 0 (7, Z) =a a g(r A 


Or OzI 
j_i-l i-1-— ; +15 
J\ (i-1\ 0 7 7 (j-2) oO" On,0(T, 2) 
9 Satay [02 817, 2)) —S 
4—1— 2 U 
=e ‘er n Fe n OT0z 
g—1 4 : 2 = a “ 
j-1l\/i\ oon oj-i-l ants, 9 TB 
a », : ( I ) (") Or?-7 = K2(7, 2) — 
=0 n= 
j-li-1.. F i j = 
eo 1 i-l gi-1—-n gj-i-l ant! é, o(r, Zz) 
. 1=0 n=0 ( l ) ( 7 ) Ori-1—-7 | A@zi-1—! nara) Or"dz! : 


(14) 


Third, we approximate w by Up,m in the rectangles Dn.m,n = 0,1,...,N-1 
and m=1,...,M —1 such that 


p-1 pitas 
1 ors nim nyom 7 } 
Un,m(T; z) = » 4! ae i - Te T™)'(z Zin)" ; (7, Z) € Drm, 


4 qj! Or Oz 
i+j=0 
(15) 
where 0,,m is the exact solution of the integral equation: 
n-1 Te41 # 
Onm(T, Z) =9(7, z) + y/ Rat, 2) Ue sq (by8 at +f Ki(t, Z)On m(t, z)dt 
€=0 TE Tn 
Mal p2zp41 z 
+ Py i K2(T, 8) Un, p(T, 8)ds +f K2(T, 8)tn m(T, s)ds 
p=0 ” 2e Zm 
malmal prega pep 
+ S- > a / K3(t, 8)ve,p(t, s)dsdt 
€=0 p=0 "75 eh 
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n-1 Te41 Zz 
+ >| / K3(t, 8)vem(t, s)dsdt 
€=0 zm 


TE 


m1 


T etl 
+ > | ‘| K3(t, $)Un,p(t, s)dsdt 
p=0 % Tn 4 2p 


4 | | ele OO cade 


(16) 


ies ; aitie, ; 
and gr tno(Tn Am) is the exact value of aa at point (Tr, 2m). 


Differentiating equation (16) j-times with respect to z and i-times with 


respect to 7, we obtain 


O99, an (7,2) 


Or OzI a 
= 01°05? (7,2) 
jot-lypn 7; i—1—r tis 
j\ (i-1\ 0 AGT I Onm(T, 2) 

EEG) n ) grav [2 Ki(t,2)] aaa 

ase j-1 i\ oT q-tt OG ke) 
PDD 4) ta) ope (pee | prog 

l=0 7=0 

5 (FN (6-2) NT gay] tema) 
. 1=0 n=0 l n ) Ori-1-n Fs OrnOz! 


4 Convergence analysis 


We consider the space L*°(D) with the norm 
Il¥llz-0(p) =inf{C € R: |y(7,z)| <C for ae. (7,2) € D} < oo. 
The following lemmas will be used in proving the convergence of the presented 


method. 


Lemma 1 (Taylor’s theorem for functions of two independent variables). 
Let f be p times continuously differentiable on D = [a,6] x [c,d] and let 
(70, 20) € D. Then for all (7, z) € D, we have 
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p-1 pis 
1 Ord f (79 Z0) i : 
=— ? J 
f(7,2) = 2 ij! = Or*OzI Po) 20) 
i+j=0 
1 Ors f (74 21) 7 ‘ 
2 J 
- x ij! = Or*Ozd (7 — 70)'(2 — 20), 
t+j=P 
where 
T = 07+ (1-8) € [a, BI, 
21 = 0z+ (1-0) € [c,d], 


Lemma 2 (Gronwall-type inequality [21]). Let w(r7,z) and p(t, z) be non- 
negative continuous functions in = [a,b] x [c,d], and let p(7, z) be nonde- 
creasing in each of the variables in Q and satisfy the following inequality: 


T 


w(r,2) <p(t,z) + i 


a 


+ «| 7 w(t, s)dsdt, (7,2) €Q, 


w(t, z)dt + a w(T, s)ds 


where « is positive constant. Then there exists a positive constant v such 
that 
w(T, 2) < vp(r, z). 


Lemma 3 (Discrete Gronwall-type inequality [5]). Let {«;}7 9 be a given 


nonnegative sequence and let the sequence {€,,} satisfyep9 < pp and 


n-1 
En Spot D> yey, n>, 
j=0 
with po > 0. Then 
n-1 
En < po exp So ky », n>. 
j=0 


Lemma 4 (Discrete Gronwall-type inequality of two variables [23]). Let wr m 
be a given nonnegative sequence, and let b;, b2,b3 and 8 be independent of h 


and k and strictly positive. If the sequence wy, satisfies 
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n-1lm-1 
enn i em tS tng ths 3D oe +B (18) 
€=0 p=0 
for alln = 0,1,...,N,m=0,1,...,M, then 
m < Bexp(y(Nh + Mk)), (19) 


where y= $ (b: + bg + 4/ (by + by)? + 1s : 


Lemma 5. Let g, &1,K2, and K3 be p times continuously differentiable on 
their respective domains. Then, there exists a positive number a(p) such 
that for alln =0,1,...,N—-—1,m=0,...,M—-—1, andi+ 7=0,1,...,p, we 
have 


< a(p), 


| OF On sm 
L&(Drym) 


Or? OzI 


where 60,0(T, z) = w(T, z) for (7,2) € Doo. 


Proof. Let oy = _ | Be IL-0 


n,m)" 


First, we have for alli+7=0,1,...,p 


iy Oy 
oar < max —— 
OT OzI 


z 


6420) cph = op) (20) 
L&°(Do,o) 


Second, for i+ 7 =0, from (13), we have for alln =1,...,N-—1 


ait sete  f 3S arn aia 


a+B=0 Tr 


T€+1 
+ C3, 2 3 ai. /-> agi” dsdt + C3, 2 7. [ ayy dsdt, 
a+B=0 


and fori+j =1,...,p, from (14), we have for all n = 1,...,N —1, 


TE4+1 


z 
Gi, at + €2,2 i. att ds 
) 


j i j-1i-1 

1 +1 41 

89 Scie 3 ott toad ; ay 0 +¢3,3 5 y ano: 
l=0 n=0 l=0 n=0 l=0 7n=0 


where, the constants c, €1,1, C1,2, C2,1, €2,2, €3,2, and c3,.3 are positive and inde- 
pendent of N and M. 


Hence, for alli+7=0,1,...,p 
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n—-1 p-l j 
aitd a+ I A 
if cota, So afi! ta.2 DT ot teaah ooh 
€=0 a+6=0 1=0 n=0 
gol. 4 a 
l 0 
hag at pale 
1=0 n=0 n=0 
n—-1 p-l j-1li-1 
+ ¢3,2hk S~ x apy +e33 5.) aly +cs2hkayy, 
€=0 a+ 6=0 1=0 n=0 
which implies that 
n—1 p-1 tga 1. 
i+j l I 
ana Setah) | Dd) agate Dy) ano, (21) 


€=0 n+l=0 n+l=0 


where Cy = C11 + c3,2k, C2 = C12 + ch + C21 + C2,9k + C33 + c3,ahk. 


Now, we consider the sequence [,, = max{a;,'g,i +j =0,...,p},n = 
0,1,...,N —1, from (21), we obtain 
n—1 itjg-1 
ae <c+cp ee ae 
Ati=0 
Using Lemma 3 with the following notations 
a n—-1 
cup = 0,3) no = c+ aph>- Vg,  ky+t = C2, 
&=0 
we obtain 
n-1 i9=1 
ti << | cteap ay Fe exp S- C2 
€=0 n+l=0 
n—-1 
< cexp(p*c2) + cp? exp(p?c2) h Se Te, 
a f20 
C3 C4 
which implies that 
n—-1 
DT, <eg3 + cah S > Te. 
&=0 


Again, using Lemma 3, for alln = 0,1,...,N—landi+j=0,...,p 


ang <Tn < egexp(aca) - 
— (22) 


a2 (p) 
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Third, for i+ 7 = 0, from (16), we have for all n = 0,1,...,N—1, m = 


1 ncughd’=<t, 
n-1 p-l m—-1 p-l 
nai m Set C1 why? S- ae, + 1haptn +c22k > S- ant? + a, kann 
€=0 a+6=0 p=0 a+f8=0 
n—-lm-1 p-1 n—-1 p-l 
+¢3,2hk S~ S° y ant? + c3.2hk \~ S- Laas 
£=0 p=0 a+B=0 €=0 a+B=0 
m—-1 p-l 
+ ¢3 ahk S » ant? + cg, ahkay*, 
p=0 a+f8=0 


and fori+j =1,...,p, we have from (17) that 


: : j-lé-t 
a+, +1 +L +1 
att! <e-taaS 3 al eS y arth +c3,3 5 y Oa 

1=0 n=0 l=0 7=0 l=0 7=0 

Hence, for alli+ 7 =0,1,...,p, 

n-1 p-l i4+j-1 itj—1 
i+ S- S~ ntl 41 +1 
o, <e+e,1h Gem + €1,2 ; Glin + Crh ; OF 
€=0 7+1=0 n+l=0 n+l=0 
m-1 p-1 i+j—-1 itj—-1 
l l l 
+ c2,2k y y ae +21 y a + C2,2k y ee 
p=0 n+l=0 7+l=0 n+l=0 
n—-lm-1 p-1 n—-1 p-l 
l l 
+ ¢3 2hk s s s ae +c3 ohk y y an, 
€=0 p=0 n+1=0 €=0 n+l=0 
m—-1 p-l 
5 5 ntl 
ab c3,2hk an p 
p=0 n+l=0 
i+j—-1 i+j—1 
41 41 
+ ¢3.3 s Oe on + c3,2hk ) Oe (23) 
7+l=0 n+l=0 


Consider, the sequence [ym = max{a?'?,,i+j =0,...,p},n=0,1,...,N— 
1;m=0,1,...,—1. Then by (23), we have 


n—1 m—-1 
ait <c+hbia > Temthbi2 >> Tap 
€=0 p=0 
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n-1lm-1 t+j—1 
+hkbisS> SY oTeptbia D. attt, (24) 
€=0 p=0 n+l=0 


where 61,1 = (c1,1 + ¢3,2k)p”, 61,2 = (€2,2 + ¢3,2h)p”, b1,3 = C3,2p" 
and b14 = C1,2 + c1ih + C21 + C2,2k + C3,3 + c3,2hk. 


Using Lemma 3 with the following notations 


ee a) 
Ei4g = an ms 


m-1 n-1lm-1 
pom eb Maa So Pe + Bho yy Dap + hkbis >) > Teo 
€=0 p=0 
kn4t = 1,4, 
we obtain from (24) that 
n—-1 m—1 
ait), <cexp(pbi 4) +h b1,1 exp(pbia) 55 Tem + kb1,2 exp(pbi4) ¥> Tn, 
Ya — eG —_—_, 7=0 
ba bi be 
n—-1m-1 
+hkby,3exp(pbia) ¥) SY > ep. 
a ie 0 pao 


It follows that, for alln =0,1,...,N—1;m=0,1,...,M—1, 


n-1 m—-1 n—-1lm-1 
Drm <ba + hbi 5 Vem + kb J) Tne + hkbs ¥) >> Teo. 
€=0 p=0 €=0 p=0 


Using Lemma 4, we obtain, for all n = 0,1,...,N —1,m=0,1,...,N—-1, 
that 
Trym < bs exp(yi(Nh + Mk)) < bs exp(yi(a + d)), 
OS (25) 
a3(p) 


where = 4 (b1 + ba + / (by + be)? + 1s : 


Hence, from (20), (22), and (25), by setting 


a(p) = max {a1(p),a2(p), a3(p)}, 


the proof of Lemma 5 is completed. 
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Theorem 1. Let g, «1, K2, and K3 be p times continuously differentiable 
on their respective ae Then equations (10), (12), (15) define a unique 
approximation uv € - 4 Dye um), and the resulting error function e(7, z) = 


w(t, z) — u(7, z) satisfies 
llellb~(p) < C(h+k)?, 


where C is a finite constant independent of h and k. 


Proof. Define the error e(7, z) on Dnjm by €n,m(T, 2) = w(T, Z) — Unym(T, Z) 
for aln =0,1,...,N—landm=0,1,...,.M-1. 


The proof is split into three steps. 


Claim 1. There exists a constant C; independent of h and k such that 
lleo,ollz(Do0) < Ca(h + )?. 


Let (7, z) € Do,o. By using Lemma 1, we obtain from (10) that 


OI nik 
Joona aa laags 
Pe 
Hence, by Lemma 5, we have 
a 
leoo(r.2)1 a(n) S> ented = 27) + ep. (26) 
eae P: 
Ya 
C1 


Claim 2. There exists a constant C2 independent of h and k such that 
llen,ollt-©(Dp,.0) es Co(h + k), 


for alln =1,...,N —1. Let (7, z) € Dn,o, we have from (13) that 
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w(T, Zz) — tn,o(T, Z) =f K2(T, 8)(w(T, 8) — tn,o(T, 8))ds 


Te44 
- >| K1(z, tee o(t, z)dt 


€=0" TE 


+ [ Ki(Z, t) (w(t, Zz) = On,o(t, z))dt 


Tn 


mal prepa pz 
- ae | K3(t, 5)e¢ o(t, s)dsdt 
£=0 TE 0 


ra | | ” it Na) — Ball aasd 


Hence, 
n—1 n—-1 
lo(r, 2) — bn,0(7, 2)| S D> hallee.oll z= (e.0) + >, PkaMlee oll n= (D¢.0) 
€=0 €=0 


+ xf \w(t, z) = On,o(t, z)|dt 
+ nf |w(T, $) — Op ,9(7, s)|ds 
0 
ae nf / lw(t, 8) _ On,o(t, s)|dsdt, 
Tm JO 


where & = max{||Kj||L©(p),4 = 1,2,3}. Then by Lemma 2, we have 


n—1 n—1 
lu(7, 2) — &n,0(7, 2) < | Do Aaillee oll n= (De0) + >, Phe |leeollz~(Deo) | Y 
€=0 €=0 


n-1 
= > hel + b)v \leg,all t= (De,0)> 
€=0 M 


which implies, by using Lemma 1, that 


llen,ollz-(Dn.o) S lle — Sn,oll + |lGn,0 — Un,oll 
n-1 side 
1 OI bn 0 44 
, j 
S Do MAillecolle~ (Deo) 4 » ag ortaer 
€=0 i+j=p 


Hence, by Lemma 5, we obtain 
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n—1 


ve ) 
llen,ollz-(Dy.o) S S hAx|leg,oll z-° (Deo) + p(k + k)?. 
£=0 


Then, by Lemma 3, we have 
a(p) . 
p! 
— ——S 
C2 


llen,ollt©(Dn.o) 


Claim 3. There exists a constant C’3 independent of h and k& such that 
llen,ml| °° (Damm) < C3(h =F k)?, 


for alln = 0,1,...,N—1 and m =1,...,M —1. Let (7,z) © Dnm. Then 
from (16) we have 


n-1 m—-1 


|w(7, 2) — Bnm(7 2) SO halleemll + D5 kellenoll 
€=0 p=0 


n—-1lm-1 


+ 37 ST hkellee.pll + a hkelle¢ mall + by hkrlen,pl 


€=0 p=0 


+ xf |w(t, z) — On m(t, z)|dt 


n 


ref lw(7, 8) — Onym(T,8)|ds 


+e ft fw w(t, 8) — On m(t, s)|dsdt. 


Then by Lemma 2, 


n-1 m—-1 
o(7, 2) — Onm(72)| SDS RL + KY [lee ml] + D5 kK(L + A)Y llenoll 
€=0 —— o=0 oS 
2 3 
n—-1lm-1 
+2 hk Hy IMleeoll 
€=0 p=0 


which implies, by using Lemma 1, that 
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llen,m|| Lt (Dy,0) <||w — Bn ym =F lon, m — Un nll 


n—1lm—-1 


< > halle mall + > kXsllenll + 5. S_ bkAallee,pl 


€=0 p=0 
1 arg be 
pe a 
i Gt! || OrtOzI 
t+J=P 
Hence, by Lemma 5, we obtain 
n—-1m-1 
llenmll| S is hrg|lee ml + > kAgllenoll + >, >, hkAallee,oll 
€=0 p=0 (27) 
2a cay 
Dp: 
Using Lemma 4, we obtain 
a(p 
lenml| < WP expla + Z)) (+ WP, (28) 
ee 


C3 


where 73 = S Ag +A3 + (r2 + A3)? +43). 
2 


Thus, the proof is completed by taking C = max{C), C2, C3}. 


5 Numerical examples 


The method presented in this paper is used to find numerical solutions to 
some illustrative examples. Our results are compared with the exact solu- 
tions by calculating the absolute error function |énm| = |w — Unjm| for all 
n= 0,1,...,N —1 and m =0,1,...,M—1, where w and v are the exact 
and approximate solution, respectively. The values of errors are computed 
for different values of p, N,M and collected in Tables 1, 3, and 4, which are 
displayed in Figure 1 for Example 1. In Table 5, the presented method is 
compared with the numerical results obtained by using the Chelyshkov poly- 
nomials method (2D-CPs) [20] and the two-dimensional block-pulse functions 
method (2D-BPFs) [22]. Moreover, the exact and approximate solution over 


the region ([0,1] x [0,1]) are displayed in Figure 2 for Example 2. The re- 
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sults in these examples confirm the theoretical estimates and suggest that 


the experimental order of convergence (EOC) is p (see Table 2). 


Example 1. Considering the Goursat problem, which is linear and homoge- 


neous in hyperbolic PDE [9], 


Aw(t,s)  Au(t, s) a Ow (t, 5) 


3 OsOt —Ot Os 


+w(t,s), (t,s) € [0,1] x [0,1], 


with initial conditions w(t,0) = e’ and w(0,s) = e*%. This equation is equiv- 


alent to the linear two-dimensional VIE defined as follows: 


Wil rR 


w(T, 2) =i(e" +e*)— 


+5 (fowaas Pornsase f° [wt s)dsdt) _ (29) 


The exact solution is w(7, z) = e7™T*. 
The numerical results for p = 3,4 and h = k = 0.05, 0.025 of the Taylor 


collocation method are presented in Table 1 and Figure 1. 


Table 1: Comparison between the approximate and the exact solution for Example 1 


(7, z) N=M=20.p=3 N=M=20p=4 N=M=40,p=3 
0.1, 0.1 1.62e — 05 1.57e — 05 4.05e — 06 
(0.2, 0.2 3.82e — 05 3.65e — 05 9.47e — 06 
0.3, 0.3 6.67e — 05 6.37e — 05 1.66e — 05 
(0.4, 0.4 1.05e — 04 9.89e — 05 2.59e — 05 
0.5, 0.5 1.56e — 04 1.44e — 04 3.80e — 05 
(0.6, 0.6 2.21e — 04 2.01le — 04 5.35e — 05 
0.7,0.7 3.04e — 04 2.74e — 04 7.34e — 05 
(0.8, 0.8 4.1le — 04 3.66 — 04 9.86e — 05 
0.9, 0.9 5.48e — 04 4.82e — 04 1.30e — 04 
(1.0, 1.0 1.99e — 03 6.39e — 04 3.38e — 04 
CPU time/sec 30.81 47.25 742.68 


Example 2. Consider the linear nonhomogeneous Goursat problem [8] 


O7w(t, 8) 
Osdt 


= dts —t?s? +w(t,s), (t,s) € [0,1] x [0,1], 
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0.000100 0.000100 
0.00009} 0.00009s} 
0.000090] 0.000090} 
0.000085 0.000085 
0.000080 ‘0.000080 
0.000075 \ by 0.000075 
0.000070 | pe 0.000704 
0.00006] 


0.000060 0.000060 
0.000055 0.00005 
0.000050 0,0000504 
0.000045 0.000045 
0.000040 0.000040 
0.000035 0.000035 
10030 
0.000025 0.000025 
0.000204 0,0000204 
0.000015 1001: 
0.000010 0.0000 104 
0.000005 0.000005 


Figure 1: Plot of the absolute error function for Example 1; left: h = k = 0.05, right:h = 
k = 0.025. 


Table 2: EOC for Example 1 


(N,M) (2,2) (4,4) (8,8) (16,16) (32,32) (64,64) 
p=2 / 1.45 1.71 1.85 1.94 1.96 
/ 2.40 2.69 2.84 2.92 2.96 


with initial conditions w(t,0) = e and w(0,s) = e*. This equation is equiv- 


alent to the linear two-dimensional VIE defined as follows: 
7323 Tf? 
w(T,z) =e" +e* +7727 — =i +f / w(t, s)dsdt. 
o Jo 


The exact solution is w(7,z) = 7227 + e™*?. 
The numerical results for p = 3 and h = k = 0.05,0.025 of the Taylor 


collocation method are presented in Table 3 and Figure 2. 


Example 3. Consider the PDE with variable coefficients 


O7w(t, 8) 9, Ow(t, 5) > Ow(t, 5) 
pipe ge oe 
+ tsw(t,s) + f(t, s), (t, s) € [0,1] x [0, 1], 


with initial conditions w(t,0) = cos(t) +e’ and w(0, s) = cos(s)+1+s7. This 


equation is equivalent to the linear two-dimensional VIE: 
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Table 3: Comparison between the approximate and the exact solution for Example 2 


(7, z) N=M=20yp=3 N=M=40,p=3 
(0.1, 0.1) 3.64e — 07 4.80e — 08 
(0.2, 0.2) L.77e — 06 2.29e — 07 
(0.3, 0.3) 4.77e — 06 6.13e — 07 
(0.4, 0.4) 1.00¢ — 05 1.28¢ — 06 
(0.5, 0.5) 1.83e — 05 2.32e — 06 
(0.6, 0.6) 3.07e — 05 3.98e — 06 
(0.7, 0.7) 4.85e — 05 6.21e — 06 
(0.8, 0.8) 7.32e — 05 9.36e — 06 
(0.9, 0.9) 1.06e — 04 1.36e — 05 
(1.0, 1.0) 1.76e — 03 2.28 — 04 

CPU time/sec 27.71 701.59 


(a) Exact solution (b) Approximate solution 


Figure 2: Numerical results of Example 2. 


w(T, z) =g(T, z) + fe + z)w(t, z)dt + fe + s”)w(r, s)ds 


+f [2 tute save 


where g(r, z) is chosen so that the exact solution is w(7,z) = cos(7 + z) + 


e7 + 27. 
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The numerical results for p = 3 and h = k = 0.05,0.025 of the Taylor 


collocation method are presented in Table 4. 


Table 4: Comparison between the approximate and the exact solution for Example 3 


(r, 2) N=M=20p=3 N=M=40,p=3 
(0.1, 0.1) 1.32e — 05 1.66e — 06 
(0.2, 0.2) 2.84e — 05 3.73e — 06 
(0.3, 0.3) 4.72e — 05 6.54e — 06 
(0.4, 0.4) 7.13e — 05 1.05e — 05 
(0.5, 0.5) 1.03e — 04 1.62e — 05 
(0.6, 0.6) 1.46e — 04 2.44e — 05 
(0.7, 0.7) 2.07e — 04 3.64e — 05 
(0.8, 0.8) 2.98e — 04 5.47e — 05 
(0.9, 0.9) 4.44e — 04 8.47e — 05 
(1.0, 1.0) 1.04e — 03 1.89e — 04 

CPU time/sec 62.96 1207.40 


Example 4. Consider the linear nonhomogeneous Goursat problem 


dPu(t,s) __ Ou(t,s) _ u(t, s) 
OsOt ot Os u(t,s)+ f(t,s), (,8) € [0,1] x [0,1], 


with initial conditions u(t,0) = e’ and u(0,s) = e?%. This equation is equiv- 
alent to the nonlinear two-dimensional VIE of the first kind [22]: 


1 
9 


T eZ 
(ge et tte et t% err) =|} | 2e7**u3 (Et, s)dsdt. 
0 JO 


It is also equivalent to the following linear 2D-VIE of the second kind: 


u(r.2)=ar.2)~ [u(t de~ fo ursyas— f° [ute s)dsdt, 


where u = w® and g(r7,z) is chosen such that the exact solution is u(r, z) = 
eTt2z, 

In Table 5, the numerical results for p = 3 and N = M = 64 of the 
Taylor collocation method are compared with the numerical results obtained 
by using 2D-CPs [20] and 2D-BPFs [22]. 
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Table 5: Comparison of the absolute errors of Example 4 


(2-*,2-*) | 2D-BPFs [22] | 2D-CPs [20] | Present method 
i=1 1.0e—1 3.5e — 5 6.le—6 
7=2 4.6e — 2 2.0e — 6 2.6e — 6 
i=3 2.9e — 2 1.5e—5 1.3e —6 
t=A4 2.3e — 2 1.2e—5 7.2e—7 
a=5 2.0e — 2 5.9e —5 3.7e — 7 
i=6 3.le—2 9.6e — 5 1.9e—7 


6 Conclusion 


In this paper, a collocation approach was presented that involves utilizing 
Taylor polynomials to find the solution to a two-dimensional linear VIE of 
the second kind, which is a conversion of a hyperbolic linear PDE Goursat 
problem. The method’s convergence and error were investigated, and several 
numerical examples were provided to demonstrate its efficiency and accuracy. 
The results showed that the method is convergent with a high level of preci- 
sion, and the numerical results match the theoretical estimates. It is recog- 
nized that this approach can be expanded and used to solve three-dimensional 


Goursat problems in linear hyperbolic equations of the third-order 


mt 3 bape = og el 


Ugjali<g 


where ¥i,5,Wi,,2 <j,1,j =1,2,3,y and F are given real functions. 
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